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' We study the cosmology of models with four space and one time dimension 

where our universe is a 3-brane and report a few results which extend existing 
work in several directions. Assuming a stable fifth dimension, we obtain a 
solution for the metric, which does not depend on any arbitrary parameters. 



We discuss some implications of this result. 



I. INTRODUCTION 
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Over the last two years there has been a lot of interest in models where our universe is 
a 3-brane (a hyper-surface) embedded in a higher dimensional bulk. The standard model 
particles are confined to the brane whereas gravity propagates in the bulk. Such models were 
conjectured early on fl|] as interesting possibilities and have recently been argued as plausible 
solutions of type I string theories ||. One of the attractive features of these models is the 
intriguing possibility that the fundamental scale, M, identified as the string scale, could be 
lower than the Planck scale, M Pe = (8tiG n)~ 1 ^ \ by several orders of magnitude ||, perhaps 
even of TeV range ||. This last possibility may provide a new way to solve the hierarchy 
problem between the electroweak scale and the scale of gravity i.e. both scales being of 
same order, there is no hierarchy to worry about. The large value of the Planck scale in this 
picture owes its origin to the existence of very large hidden extra dimensions in nature. The 
price one has to pay is of course that now one has to understand why the extra dimension(s) 
is(are) so large. The relation between the fundamental scale, the Mp£ and the volume of 
the external space, V n , is given by M 
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M 2+n V n = Mj, e . (1.1) 

This picture leads to a modification of the inverse square law of gravity at small distances, 
r ~ V^' n and can therefore be probed experimentally. If M is in the TeV range, string 
theories become accessible to collider tests. All these make the idea phenomenologically 
quite attractive. It is therefore interesting to study the cosmology of these models. 

One of the first things that one needs to know in order to study the cosmology of these 
models is the time dependence of the metric. In addition, the metric will also have a 
dependence on the bulk coordinate y, even when the bulk is totally empty, simply because 
the presence of a brane induces a nonzero curvature. It turns out that this fact leads to a 
variety of interesting consequences for the cosmology of such models j^— [14] . In particular, 



the new bulk-brane picture seems to drastically modify the standard time evolution law 
in the brane-confined universe. The five dimensional Einstein equations first studied in 
|| implies that the Hubble parameter H is proportional to the density on the brane, p, 
instead of the usual H ~ yfp of the standard big bang cosmology. Since the successes of the 
standard cosmology such as nucleosynthesis and the common understanding of subsequent 
evolution rely crucially on the assumption of H ~ -/p, a great deal of work has been devoted 
to understand how this standard behaviour of H can be recovered. Some ideas proposed 
to solve this problem include cancellation of bulk and brane cosmological constants |]10|| , 



consideration of a non zero thickness of the brane 111 etc. It has also been noted that 



Friedmann equation could be recovered on the basis of fine tuning of cosmological constants 



even when the extra dimension is not stable 121 and no matter in the bulk. However, in the 



process it has been found that a free parameter C, a constant of integration, appears in this 
equation. It contributes to the evolution of the Hubble parameter in the form of an effective 
radiation term, jeopardizing the cosmic scale parameters. The same arbitrary parameter 



has been used to study cosmological evolution even when the bulk radius is stable |L2 
We reexamine this in this paper. Writing the equations governing the evolution in time, 
t and the fifth coordinate, y of the various parameters defining the metric (the analogs of 
the Friedman- Robertson- Walker (FRW) scale factor), we show that it is possible to get a 
solution for the metric as a function of the bulk coordinate without any arbitrary constant, 
if the bulk radius is stable. This is the main result of this brief note. We then point out 
some of the implications of this result. 

This paper is arranged as follows: in section 2, we review Einstein equations in a diagonal 
metric, commonly used in the literature; in section 3, we look for solutions of Einstein 
equation keeping the bulk radius stable and obtain an explicit form for the metric which 
involves only the known FRW scale factor and no extra parameters. We recover the known 
behaviour for the Hubble parameter, i.e. H depends linearly on the brane energy density. 
We also show that the equation involving the brane-space like components can be derived 
from the other equations; they reduce to the acceleration equation for the FRW scale factor. 



Finally, we reconsider some examples already studied in |12| and rederive the exact solutions 
for the metric in the prescence of a cosmological constant in the bulk. In section 5, we make 
some remarks on the general nonfactorizable nature of this solutions. 
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II. BASIC FRAMEWORK 



The basic framework for our discussion is a five dimensional space-time where a flat (zero 
thickness) brane is localized at the position identified as y = along the fifth dimension. 
Although in this paper, we assume that the extra dimension is compact, this is not crucial 
for our conclusions and our results also apply when the extra dimension is noncompact. 

Since we are interested in the brane cosmology, we start by adopting the cosmological 
principle of isotropy and homogeneity in the three space dimensions of the brane. The 
presence of the brane clearly breaks the isotropy along the fifth dimension and this is reflected 
in the explicit y dependence of the metric tensor which we choose to have the following form: 

ds 2 = -n 2 (y, t)dt 2 + a 2 (y, t)-f ij dx i dx j + b 2 (y, t)dy 2 . (2.1) 

Here 7^ = f(r)5ij, with / _1 (r) = 1 — kr 2 being the usual Robertson- Walker curvature term, 
where k = —1,0,1; t and x 1 , i — 1,2,3 are the time and space-like coordinates along the 
brane respectively. 

The five dimensional Einstein equations take the form 

Gab = Rab ~ \g AB R = 4 [f AB + 5% 5 B 5(by) ] , (2.2) 

where k\ = 8ttG^) = M~ 3 is the five dimensional coupling constant of gravity, Rab is the five 
dimensional Ricci tensor and R the scalar curvature, A, B = 0, 1, 2, 3, 4 and /i, v = 0, 1, 2, 3. 
In conformity with the usual practice, we identify the mass parameter M with the string 
scale. In the last expression the various source terms have been explicitly separated. T^ v 
and Tab represent the stress-energy-momentum tensors of the brane and bulk respectively. 
For the scenarios we are going to discuss from now on, it is sufficient to work in the perfect 
fluid approximation to those tensors 

f A B = diag{-p B , P B , P B , P B , P T ), 

T^ u = diag(-p b ,p b ,pb,p b ). (2.3) 

where p B , P B and Pt represent the densities and pressure on the bulk and, respectively, p b 
and pb are those on the brane. Notice that assuming T 4 = avoids the complication of a 
matter flux along the fifth dimension. 

So far, most of the attention paid to this model in literature appears to have focussed 
on the case where vacuum energy (a cosmological constant) is the only component of the 
bulk stress tensor. In general however, all the components of T B could be functions of time, 
if they are in the brane and could depend on both t and y if they describe the bulk fllTf . 
Therefore, in our analysis, we will keep the energy-momentum tensor dependent on t (and 
y for the bulk terms). It is of course straightforward to see that more branes could be 
considered by including their corresponding stress tensors in Eq. ( |2.2|) and our discussion 
below easily generalises to this case. 

In order to solve the Einstein equations on the presence of the delta-function type den- 
sities, we proceed as follows. First we observe that the brane divides the bulk into two 
different domains, where the only source is Ta B - We then solve the equations in each do- 
main separately, and impose the boundary conditions at the brane to get the global solution. 
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This also helps to define the metric on the brane itself. First, the metric tensor, g, clearly 
should be continuous, i.e. the solutions must satisfy 



9ab(v = ) = g AB (y = + ). 



(2.4) 



Next, as has already been noted, since Gab involves up to second derivatives on the metric 
tensor with respect of y, we must use them to match the delta function distributions [g. 
Technically speaking this means that the extrinsic curvature K AB i n the Gauss-Codacci 
formulation JT3] should be discontinuous at the position of the branes. Then, by integrating 
Eq. (|2.2|) at both sides of the brane, we will get matching conditions for the first derivatives. 
This leads to the constraint 



o+ 



dy b G^ = k\T^ v . 



(2.5) 



To evaluate this integral we should assume that all other terms not involving second deriva- 
tives on y are finite. 

Typically, a parity symmetry P : y — > — y is assumed as in the Horava-Witten model [0 . 
Physically, this symmetry could be seen as a residual effect of the broken isotropy along 
the fifth dimension, as in the S l / Z2 orbifold construction of ref. 0. We will assume this 
hereafter. In the presence of other branes that explicitly break this symmetry, our discussions 
have to be reconsidered depending on the number of branes. 

Let us now proceed to the details. Using the above form of the bulk stress tensor, the 
non trivial components of the Einstein equations (away from the brane) are given as || 



G 



00 



Gij 



Gq4 
G44 




K 2 5 n 2 p B , 



Hi + 

lij - = Kla 2 P B jij, 



4b 2 P T ; 



(2.6) 



(2.7) 
(2.8) 
(2.9) 



where primes (dots) are used to denote derivatives with respect to y (t). This system of 
equations is supplemented by the Bianchi identity T A B . A = 0, which translates into the 
conservation laws [TT| 



PB + 3- (p B + P B ) J r h T (PB + P T ) 

a b 

o! Til 

P > + 3 _ (p T -P B ) + - (P T + p B ) = 
a n 



(2.10) 
(2.11) 



Next, by using Eq. (|2.5|), the following boundary conditions are easily obtained 
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A "' (2.12) 



ab 
An' 



nb 



2 
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3 (3p 6 + 2p b ); (2.13) 



where the left hand side of the above equations has to be evaluated at the position of the 
brane, and the function Aa'(O) := a'(0 + ) — a'(0~) give the size of the jump of the derivative 
of a(y). The same applies for An'. Since we are assuming P parity, the jump on the above 
equations could be expressed in terms of the limiting value on one side of the brane, for 
instance by Aa'(O) = 2a'(0 + ), and a similar relation for An'(O). 

Let us notice that if we use the above boundary conditions, we may evaluate equation 



( p.8| ) on the brane to get the conservation equation 



a 



p b + 3-(p b + Pb) = 0. (2.14) 



a 



This is a general result that is independent of the bulk content. In the subsequent discussion, 
we will assume that the bulk is stable and use the freedom of coordinate redefinition to set 
6 = 1. This condition will simplify the Einstein equations making it easier to extract its 
physical meaning, as we will see later. 



III. MASTER EQUATIONS OF BRANE COSMOLOGY 



Once we assume that the fifth dimension is stable, we can follow standard procedure as 
in four dimensional cosmology to reduce the equations given in the last section to a minimal 
set. First, notice that the conservation laws (|2.10|) , (P-ll|) and (|2.14f) , will have 6 = 0. Taking 
6 = 1, we notice that the Goo component (Eq. |2.6| ) of the Einstein equations is already the 
equivalent of the Friedmann equation, with the Hubble parameter defined as a function of 
y. However, the presence of the brane will require that the term with second derivative be 
regularized by extracting the divergent part. The Friedmann equation then has the form 



H 2 (y) := - 



= pB + n 
a J 3 




n 2 
a z 



(3.1) 



where a R stands for the regular part of the function. Clearly, we recognize an expression 
similar to that given early [12]. However, let us stress that there is no unknown constant of 



integration as in [|12|]. Also, we may identify the regular term as the contribution of the Weyl 
tensor of the bulk |13|]. Since this expression is continuous, thanks to parity symmetry, we 
may evaluate it on the brane to get the effective Friedmann equation of our universe, where 
the time component of the metric tensor is chosen to be no = 1: 



Hi = ^Pbo + ( yPfe 




k 

2 ' 



(3.2) 



where the subindex stand for the evaluation at y = 0. As already known ||, this expres- 
sion has the squared dependence on the brane density. On the other hand, it also has a 
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dependence on the metric outside (i.e. in the bulk), through a B , which however could be 
evaluated, once we solve for a as a function of y. Particularly, for simple cases as an empty 
bulk, the y dependence can be explicitly extracted and this term evaluated. 

Next, let us consider the equation involving Gij, (p/7| ) . Again, we separate the singular 



and the regular parts of the second derivative term and by introducing ( |3.1| ) we reduce this 
equation to a form equivalent to the acceleration equation given by 



a o 




+ 



n 



R 

n 



(3.3) 



Note that, as in the standard FRW cosmology, this is not an independent equation. Indeed, 
it can be derived by taking the time derivative of ( |3.1|) and combining that expression with 
the energy conservation law ( |2.1(JD and G04 equation (|2.8|). This derivation holds regardless 
of whether the bulk radius is constant or changing with time. As a result, this equation 
does not provide any extra information, but it will be useful in what follows. 

We now turn to the equation involving G 44 , ( |2.9| ). In conjunction with the equation 
involving Goa, this represents the new ingredient of the brane cosmology and can be a 
window to understand the y dependence of the cosmological parameters. Notice that our 
procedure is in contrast to that used in previous works where the Friedmann equation has 
been obtained from G 44 component. If we substitute Eqs. ( |3.1| ) and ( |3.3[ ) in the equation 
for the G 44 component (Eq. (2.9)), we find a simple equation that governs the behaviour of 
a(y,t) and n(y,t) on the bulk (into each one of the domains). It is given by 



3^ + ^- 

a n 



Hr. 



(3P B - 2P T - p B ) 



This expression is supplemented by G 4 which trivialy reduces to 

n(y,t) = \(t)a(y,t), 



(3.4) 



(3.5) 



where X(t) is an arbitrary function of time. Using the freedom of fixing the gauge on the 
coordinate system, we can set no = 1 in which case we get A = do ■ However in most of the 
results this choice is not necessary at all. 

Let us emphasize that Eqs. ( |3.1| ), (|3.4|) and ( |3.5| ) form the set of master equations, in 
the sense that they determine the t and y dependence of the metric. We study the solutions 
of these equations in the next section. 



IV. EXACT SOLUTIONS ON THE BULK 

As a simple application of our master equations let us reconsider the cases already studied 
in the literature. Let us assume that the stress tensor of the bulk gets contribution only 
from a cosmological constant, T A B = —A B 5 B . Then, Eq. (|3.4| ) reduces into 



a" n" 2 „, , , 

3- + - = --4K B . 4.1 
an 3 



This equation, together with the scaling equation fl3.5|) , can be solved in a straightforward 
manner. For A B = we get a linear solution just as in || 
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a(y,t)=A\y\ + B; n = X (A\y\ + B) . (4.2) 

Here we have already imposed P (parity) symmetry on the solution. Next, by using the 
boundary conditions (|2.12Q and fl2.13|) we get the final result 

a(y,t) = ao (l-y^bl); and n{y,t)=no U + ^(3p b + 2p b )\y\). (4.3) 



Notice that in the expression for n requiring consistency with the scaling ( |3.5| ) leads to the 
conservation law 

p b + 3H {p b + p b ) = 0. (4.4) 

Note that a'^ = in this case, and there is no obvious way to get the correct Friedmann 
equation ( |3.2| ) (i.e. linear rather than squared dependence of H 2 on p b ). 

Next, we assume that A# is non zero. Again ( |4.1|) is easy to solve, and after using the 
boundary conditions, we get for A^ < 

a(y,t) = a (^cosh(p\y\) - ^-p b sinh(/x|y|)^ , 

n(y, t) = n (cosh(p\y\) + ^{3pb + 2p b ) smh(p\y\)\ ; (4.5) 



where p? = —kIAb/6. We recognize those solution presented by Binetruy et al. in fi~2| 



when they take their integration constant as zero. Our finding is that this is the only 
allowed solution with a stable extra dimension. In this case, the condition ( |3.5| ) reduces to 
( |4.4|) . For Ab > the solutions have a similar form, with the hyperbolic functions replaced 
by cos and sin respectively. 

If we also assume a brane tension, A^, the Friedmann equation on the brane becomes 

K t A „ , K t „2 k , K 5 A 2 , K 5 A ^ q\ 



m = -^A bPb + -±pi n ~\~ -^Ai + —As- 



We emphasize that this result shows that C, the unknown constant found in |T2[ is actually 
zero. A fine tunning among the cosmological constants [10] 

A* = -f A? (4-7) 

may linearize the Friedmann equation in the limit where p b A&. As it is clear, this will 
work only if the vacuum energy of the bulk is negative. Moreover, to get the right expression 



of the Friedmann equation as in standard cosmology, an extra fine tunning [10 



A b = 6^ (4.8) 



is required, where n\ = M p \. With those new assumptions the solution ( |4.5|) reduces into 
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a = a (e-^l-^sinh^lyl)) 



n = n (e~ M + 3pb + 2pb sinh(^|y|)) (4.9) 

Clearly, if we neglect the contribution of the brane densities, we identify the Randall- 
Sundrum static solutions 



a = aoe 



' M and n = n e- M . (4.10) 



So far we have only been concerned with finding solutions of the five dimensional cosmo- 
logical equations with an empty bulk. If we wanted to include a bulk field, e.g. to provide 
a new picture of inflation |15| , one would have a modified form of the Equation |3.4j: 



n" n" k 2 

3- + - = -2^(V(0)+A B ). (4.11) 
an 6 

Clearly, if the scalar field evolves very slowly, e.g in the inflation epoch, the solution will be 
approximately of the form ( f4.5| ), but now with 

^ = ~4 (V(<f>)+A B ). (4.12) 
o 

By setting the fine tunning ( f4.7| ) on ( |3.2j ) we get 

Hl ^l v[m) ^yA^A pl (4 . 13) 



As expected, the inflaton potential contributes linearly to last equation, and the effective 
potential is just the value of the bulk potential on the brane. 



V. NON FACTORIZATION OF THE METRIC PARAMETERS 

Before closing the present discussion, let us comment briefly on the nature of the exact 
solutions. First we point out that all the solutions presented on the previous section are not 
factorizable. This is in fact a general property of this class of models as we will show now. 
For this purpose, let us start with case when the bulk radius b is time dependent and we 
will see that factorization of the metric parameters n(y,t) and a(y,t) will not be consistent 
with a stable b. For this purpose, we start with the identity 




a' a a' 
a a a' 



(5.1) 



which, together with the Gq4 equation leads to 






a 



(5.2) 



Note that this expression reduces to Eq. (|2.14|) on the brane. 



S 



Now, let us assume that a and n are factorizable i.e. 



a = ao(t)pi(y) n = n (t)f3 2 (y). (5.3) 
One can then easily see that 




(5.4) 



This can be used to rewrite Eq. (|5.2|) as 




P1P2 



1=0. (5.5) 



This equation can be integrated in a straightforward manner and leads to 

u _ „-&(»). 



(5.6) 



where is given in terms of fti^ by the expression between parenthesis in ( |5.5|) . From this, 
we conclude that the factorization ansatz works when the bulk radius, b is time dependent 
since clearly a grows with time. A stable b would then mean that f3% must equal zero, which 
is possible only if (3\ = fa- Furthermore, since Eq. ( |5.4|) holds everywhere, we may evaluate 
it on the brane and then conclude using Eq. (2.12) that p\ = 0. From this, it follows that 
factorizability of n(y, t) and a(y, t) proposed above implies that only a cosmological constant 
may be present in the brane. Since we require a true time dependent density in the brane 
to describe realistic cosmology e.g. the transition from a radiation to matter dominated 
universe, the exact solution to the metric can not be factorizable. 

Parenthetically, let us note that if we had chosen a more general form for the exact 
solutions 



a = a (t)(3i(y, t) n = n (t)/3 2 (y, t); (5.7) 

where the functions (3\^ satisfy the boundary condition /?i(0,t) = /?2(0,t) = 1, no such 
restriction on brane energy density would emerge. This is indeed the form of the solutions 
discussed above. For this case, the five dimensional scalar curvature can be written as 

R = /? 2 - 2 %) + • • • ; (5.8) 

where is the four dimensional scalar curvature formed by a and n , the dots represent 
extra terms. 

At this point, we note a puzzling feature with regard to the true definition of the New- 
ton's constant. One can define the Newton's constant in two ways: one by identifying the 
coefficient in front of the energy density in Friedmann equation and another way is by inte- 
grating over the fifth dimension in the action integral. We call the first definition a "local" 
one whereas the second one we call "global". Clearly in the presence of matter, the first 
one gives a constant Gn whereas the second gives a time-dependent Gn- We can check this 
easily as follows. Consider the gravity action in five dimensions 



9 



S = Jd 4 xdy V=9^R, (5-9) 

Using the relation 

^=^l(3l(3 2 (5.10) 

(with <?4 the determinant of the corresponding four dimensional metric) and Eq. (5.8), we 
obtain 

- 2 =\ fdy{3lfc l (5.11) 

adopting the "global" definition of Newton's constant. The integral in the last equation is 
taken over the whole fifth dimension. It is clear from the above expression, that this leads, 
in general to a time dependent Newton's constant whereas if we used Friedmann equation, 
we would get it to be time independent. We wish to note that when one uses the static 
solution, a la Randall and Sundrum, there is no time dependence due to the absence of 
matter and hence no puzzle. 



VI. CONCLUSIONS 

In summary, we have analyzed Einstein equations for cosmology in five dimensions within 
a brane-bulk picture. Restricting to the case of a stable bulk radius, we extract the general- 
ized Friedmann equation which does not contain any integration constant, but a term that 
involves a regular part of the second derivative over the spatial component of the metric 
along the brane coordinates. We then evaluate this second derivative using the complete 
set of the master equations for five dimensional cosmology and find that for the case of a 
stable bulk radius, there is no arbitrary constant in the Friedmann equation in the brane. 
This makes it easier to interpret the brane cosmology as the standard big bang picture. An 
advantage of our analysis is that we do not need to make any explicit assumptions regarding 
the bulk content. In this sense it generalizes the results presented before in the literature 
where only a cosmological constant was assumed to be present in the bulk to obtain the 
Friedmann equation. Finally we make some remarks on the general nonfactorizable nature 
of the exact solutions. 
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